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SECTION  I 


INTRODUCTION 

./ 

This  report  is  ottered  for  standardizations  in  air-to-air  gunnery  assess- 
ment modeling  (ATAGAM) . While  directed  toward  this  charter  it  is  also 
applicable  to  all  forms  of  aerial  weapon  delivery. 

The  treatise  evolved  from  an  original  intent  to  clarify  certain  aspects 
of  specific  ATAGAM.  Clarification  required  definition  of  the  basic  frame- 
work. Suggestions  from  associates  expanded  it,  so  that  only  a few  additions 
were  necessary  to  achieve  the  present  form.  It  is  the  apparent  solution  to 
the  quandary  of  documenting  each  and  every  model  which  uses  essentially  the 
same  transformations  as  about  12  other  models  within  the  IJSAF  ATAGAM 
commun i ty . 

There  are  24  different  ways  of  indexing  an  orthogonal  reference  frame 
(RF)  along  the  basic  directions  of  3-dimensional  Euclidean  space;  i.e.,  the 
six  basic  directions  for  one  axis  times  the  four  cardinal  points  of  rotation 
for  another  axis.  Directions  of  two  axes  of  an  orthogonal  RF  defines  the 
third.  For  each  indexing,  there  are  six  different  ways  of  oid>  ring  the 
rotations  of  heading/bearing/psi , elevation/pitch/theta,  and  roll/bank/phi 
(throughout,  the  Greek  notations  will  be  generally  used  to  denote  specified 
rotations  after  associations  with  common  names  have  been  made)  There  are 
also  six  different  transformations  for  the  three  primary  reference  frames 
of  flight;  i.e.,  inertial,  body,  and  wind.  When  considerations  are  made 
for  a variety  of  secondary  reference  frames  and  multiple  bodies  of  interest, 
it  is  readily  apparent  that  bookkeeping  is  a major  concern.  Standards  of 
indexing,  rotational  ordering,  and  inter-reference-frame  transformations 
that  minimize  bookkeeping  are  preferred. 

Conciseness  has  been  subordinated  to  an  attempt  at  thorough  explanations 
to  beginning  Scientific  and  Engineering  personnel.  Knowledge  of  right  hand 
and  left  hand  RFs  including  indexing  and  sensing  is  assumed.  The  great 
majority  of  mathematical  operations  in  the  development  is  defined  by  matrix/ 
linear  algebra  with  occasional  inclusions  of  vector  dot-product  and  cross- 
product  operations. 

To  simplify  report  preparations,  context  is  used  to  define  vectors  and 
matrices,  instead  of  the  bold-faced  type  convention.  Vector  components  are 
identifiable  by  the  addition  of  X,  Y,  or  Z to  the  vector  symbol.  Dot  and 
cross  product  operations  are  denoted  by  (•)  and  (x)  , respectively,  between 
the  vector  symbols.  Also,  the  sub  ( ) notation  is  used  to  save  lines  t1 a* 
would  be  lost  by  using  subscripts.  No  distinctions  have  been  made  in  the 
use  of  upper  and  lower  case  use  of  x,  y,  and  z. 


ccNiiRAi.  background 


In  analytical  studies  and  associated  modeling,  confusion  sometimes 
arises  concerning  the  definition  of  the  direction  of  linear/ angular  quan- 
tities or  of  transformations,  or  both,  that  influence  gunnery  solutions. 

In  some  cases,  the  modeler  may  take  equations  from  a text  or  report  that 
may  not  be  either  correct  or  directly  applicable  to  the  problem  at  hand. 

The  following  discussions  are  offered  to  provide  both  insight  and  a basis 
for  standardization. 

for  development,  the  references  are  considered  excellent  treatises, 
providing  good  fundamentals.  Unfortunately,  definitions  and  purposes  may 
be  somewhat  vague  to  those  with  lesser  experience.  Also,  certain  authors 
may  use  unique  definitions  which  cause  undue  complications  when  merged  with 
overall  modeling.  This  effort  is  directed  for  the  consistency  which  avoids 
these  complications. 

The  fundamentals  should  be  carefully  applied,  since  errors  are  ever 
present.  Reference  1 uses  numerous  restricted/limited  definitions  that  are 
not  correct  for  general  accounting.  Hxamples  are  included  in  Section  IV. 
The  pz  of  the  v-sub  x equation  (.5.1,  b)  of  Reference  2 should  be  qz.  Also, 
the  5 in  the  z-sub  V,  z-sub  3 axis  of  figure  4.9  in  Reference  2 should  be  2 
further,  in  the  discussion  on  stability  axes  (4.2.7),  it  may  be  interpreted 
that  the  X-axis  must  lie  in  the  aircraft  (xz)  plane  of  symmetry,  while 
Reference  1 requires  it  to  be  parallel  to  the  wind  vector.  Obviously,  an 
aircraft  can  he  stable  with  a steady-state  sideslip/beta  angle.  The  author 
restricted  definition  can  be  interpreted  as  an  omitted  step  to  direct  cor- 
relations with  vertical  (xz)  plane  stability  derivatives.  The  sideslip 
transformations  in  References  3 and  I are  erroneous. 

References  1 and  2 define  beta  as  measured  from  a fixed-body  x axis  to 
flu'  xy-plane  wind.  To  be  consistent  for  the  RF,  it  should  be  shown  wind- 
to  body,  similar  to  the  ang 1 c-of-at t aek  (AOA)/alpha  definition.  Then  it 
would  be  obvious  that  beta  is  oppositely  sensed  to  psi.  Of  course,  all 
flight  instrumentation  including  wind  vanes  is  affixed  to  the  airframe; 
and  wind  angle  components  can  he  consistently  shown  to  the'  airframe.  This 
discussion  illustrates  an  example  of  mixed  sensing. 

References  1 and  2 use  !’,  Q,  and  R for  angular  rates  about  the  local  RF 
(body  or  wind).  While  Reference  1 uses  p,  q,  and  r as  smal l -perturbs t ion 
(delta)  change  in  steady- state  values,  Reference  2 uses  p,  q,  and  r as 
absolute  rates  that  are  corrected  for  both  earth  rotation  and  navigational 
changes.  Ihe  most  commonly  used  practice  over  the  years  of  using  p,  q, 
and  r for  body/wind  local  rates  is  proffered  as  standard;  I’,  Q,  and  R being 
used  alternately  (and  obviously  for  all-capitals  Fortran  coding).  Primes 
or  other  distinctions  should  be  used  for  special  case  considerations. 


Finally,  the  RF  definitions 
author  intent;  but  there  are  so 
not  relate  to  the  general  order 


of  all  references  are  clear  and  concise  for 
many  which  address  special  cases  that  do 
of  accounting  for  aerial  gunnery  problems. 


Reference  1 presents  a most  comprehensive  description  of  aerodynamic 
forces  and  moments  up  through  the  disturbed-stability  RF  while  Reference  2 
describes  origins  of  RFs  and  inter-RF  developments  with  the  modern  matrix/ 
linear  algebra  approach.  References  3 and  4 are  masterworks  of  overview 
in  exterior  ballistics  developments,  though  neglecting  relative  geometry 
and  other  simple  codes. 


Before  attempting  to  uncover  the  confusing  ramifications  of  the  various 
RFs,  it  appears  warranted  to  discuss  the  biggest  problem  that  generally 
submerges  the  subleties  of  some  interrelations.  The  major  cause  of  con- 
fusion and  mistakes  is  the  mixing  of  Rfl  and  LH  RFs,  or  indices,  followed 
closely  by  using  different  ordering  of  rotations  from  the  standard  psi- 
theta-phi.  While  any  consistent  system  is  correct,  the  bookkeeping  is 
generally  overwhelming.  No  generality  is  lost  by  defining  each  RF  as  Rll  and 
indexing  it  back  to  its  primary  RF.  This  preferred  method  will  minimize 
the  number  of  transformations. 


The  proliferation  of  mixed/erroneous  RFs  has  generally  been  caused  by 
the  desires  of  the  modeler  to  conform  with  older  conventions  such  as  the 
U.S.  Army  Artillery  convention  of  X-East,  Y-North,  and  Z-Up  and  the  gunner’s 
plane  (impact  plane)  where  high  and  right  misses  were  defined  as  positive  Y 
and  X misses,  respectively.  These  conform  to  the  lower  mathematics  graphical 
plot  convention.  Mixed  sensings  or  conversions  were  necessary. 

In  engineering  practices,  X-North/forward,  Y-East/right , and  Z-GVD 
(gravity  vector  directed)/down  has  the  obvious  attribute  of  agreeing  with 
conventional  psi,  theta,  and  phi  definitions  and  alignment  with  the  external 
force,  gravity.  North,  East  and  GVD  are  associated  with  accounting-RF 
indices,  while  forward,  right,  and  down  are  associated  with  rotatable  RFs 
of  the  objects  of  interest. 

Interpretations  in  the  old  conventions  can  be  readily  accomplished 
after  standard  RH-transformations  if  deemed  necessary.  When  used,  clear 
notation  of  intent  is  recommended. 


It  is  noted  that  a projectile/bullet  that  flies  through  the  atmosphere 
is  an  aerodynamic  body  that  is  affected  quite  similarly  to  air-craft  with 
the  exception  that  aircraft  have  variable  controls/surfaces  and  power-plants 
to  maneuver  and  sustain  energy/flight.  Thus,  RFs  for  accounting  of  the 
bullet  states  should  be  similar  but  simpler  to  those  for  aircraft.  Invariant 
configuration,  ballistic  flight,  and  no  secondary/internal  RFs  greatly  sim- 
plify the  accounting  system.  However,  there  may  be  as  many  as  30  or  more 
accountable  (i.e.,  between  firer  and  target  fly-by)  projectiles,  requiring 
stored  state  arrays  in  revolving  tables  that  drop  those  beyond  fly-by  and  add 
newly-fired  ones.  Generally,  nominal  trajectories  are  computed  for  initi- 
ations at  the  beginning  of  each  iteration  interval  and  interpolations  are 
made  in  accordance  with  the  firing  rate. 


In  concluding  the  general  background,  it  is  noted  that  all  modeling  is 
relative.  Though  some  factors  are  known  or  suspected  to  vary  in  the  real 
world,  models  neglecting  these  factors  are  chosen  where  they  will  hopefully 
have  negligible  impact  on  the  results.  Objective  interpretations  should 
always  be  qualified. 


SECTION  III 


FIXED  AND  MOVING  INERTIAL  REFERENCE  FRAMES 

The  absolute  inertial  RF  has  been  defined  as  the  center  of  the  universe 
which  is  academic  and  undefinable.  The  galactic-inertial  RF  is  of  the  same 
stature,  but  manageable  in  a restricted  sense.  A geocentric  (earth-center 
origin)  inertial  RF  can  be  defined  for  multi-body  orbital  flight  with  the 
Z directed  toward  the  sun  and  the  xz  plane  defining  the  earth's  orbital 
plane.  For  atmospheric  flight,  some  simplified  models  consider  the  earth  as 
a fixed  spheroid,  correcting  the  solution  to  varying  degrees  of  sophistica- 
tion for  earth  rotation  and/or  gravity  vector  anomalies  afterward.  Some 
sophisticated  models  consider  rotation,  oblateness,  and  curve-fitted  mag- 
netic anomalies;  but  obviously  avoid  the  mechanics  of  rapidly  varying 
anomalies  where  platforms  are  described  to  go  wild.  The  rotating  geocentric 
RF  has  its  Z-axis  directed  through  the  South  spin  pole  and  its  X-axis  di- 
rected through  the  Greenwich  Meridian.  The  earth-surface  RF  is  directly 
applicable  to  the  navigational  problem.  Initialized  at  home-base,  solutions 
determine  aircraft  state  relative  to  the  earth  and  to  points  of  interest; 
i.e.,  destinations.  Sophisticated  internal  mechanizations  are  included  to 
direct/tune  X-North,  Y-East , and  Z-GVD  and  to  compensate  for  earth  rotation. 
These  are  important  to  space  flight  or  long-term  navigation,  but  are  ob- 
viously overly  complex  to  the  few  seconds  or  less  characteristics  of  the 
gunnery  pass. 

An  important  one  for  ATAGAM  is  the  local  attitude- indexing  RF;  i.e, 
the  true  vertical  if  and  directionality  sensed  by  either  the  Vertical- 
Directional  Gyro  (7G-DG)  Set  or  the  Inertial  Reference  Platform.  Note, 
these  are  the  true  indices,  not  the  sensor  outputs  nor  the  attitude  outputs/ 
relations  of  the  aircraft/weapon  body  RFs  to  these  indices,  discussed  in 
the  next  section.  This  is  the  moving-base  inertial  RF  (MIRF)  of  the 
references . 

The  accuracy  of  the  VG-DG  Set  with  drift  rates  of  up  to  20  deg/hr  for 
accelerated  flight  can  significantly  affect  some  long-range  gunnery  problems, 
but  most  contemporary  ones  include  local-gravity  erection  or  caging  tech- 
niques which  obviates  this  concern  for  reasonable  flight  profiles,  while 
absolute  directionality  is  generally  not  a significant  factor  in  aerial 
gunnery.  Also,  the  RF  is  usually  either  physically  located  very  close  to 
the  aircraft  center  of  gravity  (eg)  or  is  compensated  for  the  offset. 

It  is  noted  that  the  origin  of  all  attitude  RFs  for  ATAGAM  is  the  body 
eg.  Variability  of  the  eg  is  discussed  in  the  next  section. 

To  account  for  the  relative  geometry  of  attacker/firer , target,  and 
projectiles,  it  is  a common  modeling  practice  to  freeze  the  MIRF  in  inertial 
space  at  significant  events  during  the  gunnery  pass.  The  primary  one  is 
used  to  account  for  firer,  targets,  and  bullet  attitudes  and  positions  and 
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therefore  relative  geometry.  Its  origin  is  fixed  at 
pass  at  the  firer's  eg  at  the  altitudo/negative-z  of 
is  initialized  with  respect  to  the  firer  in  this  RF. 
earth  approximation  is  sufficient  for  aerial  gunnery 
fixed  inertial  RF  (FIRF)  used  for  basic  accounting, 
parallel.  Both  are  called  earth  in  some  models,  overview  being  necessary 
for  discernment. 


the  initiation  of  the 
interest;  the  target 
Obviously,  the  flat 
ranges.  This  is  the 
Note,  MIRF  and  FIRF  are 


In  the  more  sophisticated  models,  projectile  fly-out  modeling  is  incor- 
porated and  compared  with  stored  target  states  to  define  fly-by;  i.e., 
intersection  in  the  impact  plane.  The  most  commonly  used  fly-out/exterior- 
ballistics  code  is  the  point-mass  (PM)  model  which  assumes  zero  lift;  i.e., 
no  yaw  of  repose.  Relative  errors  to  the  full  six  degree  of  freedom  model 
are  approximately  in  the  few  mils  or  less  category  for  design  gunnery;  i.e., 
1000  to  1500  feet  line  of  sight  (LOS)  range  and  wind  angles  up  to  about  10 
degrees.  Because  of  the  simplification  of  the  PM  model,  a single  plane 
coordinate  system  is  sufficient.  Some  PM  models  operate  open  loop,  iterating 
projectile  travel  in  the  X-forward  and  2-down  from  drag  (coefficient)  versus 
mach  tables  and  flight  equations,  while  others  use  closed  form  equations  such 
as  the  P-Q  Siacci  RF  where  P is  forward  along  the  line  of  departure  and  Q 
is  down.  Of  course,  the  simplest  model  is  the  straight  line  with  average 
velocity  over  the  ranges  of  interest.  Average  velocity  approximations  are 
available  for  both  inertial  and  relative  ranges.  Gravity  drop  of  the  no-lift 
body  can  also  be  easily  calculated.  Considering  the  gross  approximations  of 
the  PM  model  and  inclusion  of  aiming  statistics  that  submerge  exact  details, 
there  are  branches  of  the  aerial  gunnery  problem  where  simplified  ballistics 
should  be  given  consideration.  Originated/frozen  at  muzzle-exit  or  time-of- 
fire,  these  are  the  exterior  ballistics  RFs  (EBRFs) . The  single-plane  F.BRFs 
require  transformations  to  the  FIRF  fo*-  impact-plane  determinations. 


SUCTION  IV 


AIRCRAFT  FIXED-BODY  REFERENCE  FRAMES 

The  primary  fixed-body  (FB)  RE  to  which  all  other  EBREs  should  be  re- 
lated is  the  one  that  is  used  to  reference  its  attitude.  CG-originated 
and  initialised  with  respect  to  the  MIRE,  the  forward  extension  of  the 
fuselage  reference  line  (ERL),  Y-right,  and  Z-GVD  for  level  xy  plane  define 
this  body  RE  f BRF ) whose  orientation  defines  Euler  attitudes.  The  Y-right, 
a principal  axis  defined  in  Section  V,  is  invariant  in  most  EBREs. 

Obviously,  the  eg  varies  with  fuel /armament  loads  and  configuration; 
i.e.,  different  internal  gear/crew  and/or  variable  states  such  as  tilt/swing 
surfaces,  etc.  Thus,  the  FRL  must  be  defined  for  either  a nominal  or 
critical  condition  such  as  landing.  Because  weight  and  balance  (WB)  is  so 
critical  to  flight  control,  fighter  aircraft  are  designed  such  that  the  eg 
is  constrained  to  within  a few  inches  about  nominal  in  the  fore  and  aft 
direction,  less  in  the  vertical,  and  even  less  than  that  laterally;  i.e., 
fuel  and  stores  management  are  designed  to  cause  minimum  eg  shift.  Thus, 
while  eg  shift  from  nominal  must  be  acknowledged,  its  effect  on  the  BRF  can 
be  considered  negligible  to  the  aerial  gunnery  problem. 

In  some  instances,  the  modeler  may  be  faced  with  converting  eg  and  ERL 
from  the  design/WBRF  (weight  and  balance  RE)  to  the  BRE.  The  origin  of  the 
WBRE  is  generally  located  outside  of  the  aircraft  to  keep  most  dimensions 
positive;  i.e.,  in  front  of  the  most  forward  tip/nose  and  either  up-down 
located  with  an  extension  of  an  FRL-type  line  or  the  ground  plane  for 
clearance  considerations.  X,  Y,  Z dimensions  in  inches  are  given  in  fuse- 
lage stations  positively  aft,  buttline  (inside  fuselage/empennage)  or  wing 
stations  positively  right-side  looking  forward,  and  waterlines  positively 
up.  While  buttline  wing  station  is  measured  from  the  symmetrical/xz  plane, 
both  origin  and  direction  arc  arbitrary  for  the  x/z  axis,  though  an  orthog- 
onal RF.  Most  aircraft  do  not  touch  down  or  rest  level  on  their  landing 
gear.  Thus,  specifically  when  this  RE'  is  defined  to  show  ground  clearances, 
the  slope  of  the  FRL  with  respect  to  a reference  waterline  is  required. 

Figures  1 and  2,  taken  from  Reference  1,  are  presented  to  show  the  Wing 
Chord  RE  (WCRF)  and  associated  angle  definitions.  It  is  an  important  fixed 
body  RF  (FBRF) , because  lift  occurs  along  negative  z.  (WC,  called  so  be- 
cause most  lift  is  derived  from  the  wings,  is  defined  by  the  total  lifting 
area.)  Note  that  the  pitch  attitude  of  the  WC  is  not  the  primary  pitch 
attitude  of  the  BRE'.  This  frame  is  identical  to  the  BRE  except  for  a small 
rotation  about  the  identical  y-axis.  Thrust  and  Drag  REs  are  obviously 
similarly  small  rotations  of  the  BRE'.  Usually,  drag  is  accounted  for  in 
the  stability  axes  system. 

Figure  3,  also  from  Reference  1,  but  corrected  for  the  aforementioned 
beta  definition,  shows  another  restricted  definition;  the  associated  angles 
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Plane  of  Symmetry 


LEGEND : 

U - Aircraft  Velocity  Vector 
a - Angle  of  Attack 
8 - Angle  of  Sideslip 
y - Angle  of  Climb 
0 - I’ itch  Angle 
ij;  - Heading  Angle 

Subscript  0 - Denotes  steady-state  value 

slOTE:  Vertical  Plane  Angles  are  shown  for  wing  chord  line  reference 
(Figure  2).  Pitch  is  referenced  to  the  fuselage  reference  line. 


Figure  1.  Body  Reference  System 


All  Angles  Pictured  Are  Positive 

1 


LEGEND: 

U - Aircraft  Velocity  Vector 
a - Angle  of  Attack 

6 - Angle  of  Sideslip 

Y - Angle  of  Climb 

9 - Pitch  Angle 

ij;  - Heading  Angle 
Subscript  0 - Denotes  steady-state  value 

'10TE:  Vertical  Plane  Angles  are  shown  for  wing  chord  line  reference 
(Figure  2).  Pitch  is  referenced  to  the  fuselage  reference  line. 

Figure  2.  Vertical  Plane  Geometry 


LEGEND 


U - Aircraft  Velocity  Vector 
a - Angle  of  Attack 
3 - Angle  of  Sideslip 
Y - Angle  of  Climb 
0 - Pitch  Angle 
- Heading  Angle 

Subscript  0 - Denotes  steady-state  value 


NOTE:  Vertical  Plane  Angles  are  shown  for  wing  chord  line  reference 
(Figure  2).  Pitch  is  referenced  to  the  fuselage  reference  line. 


Figure  7>.  Horizontal  Plane  Geometry 


in  the  horizontal  plane.  Psi  is  shown  to  be  sensed  relative  to  the  flight 
path.  This  is  only  true  when  the  flight  path  is  North-directed ; and,  in 
general,  the  flight  path  may  have  any  orientation.  Thus,  the  psi  in  the 
figure  should  be  more  aptly  labeled  delta-psi.  Note  that  z is  into  the 
page.  Accordingly,  g is  positive  for  left-hand  rotation. 

The  fixed-sight  RF  (FSRF)  originates  on  the  sight  combining  glass.  While 
X may  be  parallel  to  BRF  X in  certain  vehicles,  it  may  be  along  the  primary 
armament  datum  line  (ADL)  in  others  or  multiple  FSRFs  may  be  defined  for  the 
various  armaments.  The  ADL  specifies  the  nominal  line  of  departure  at  the 
design  point  altitude  and  airspeed  for  missile,  rocket,  and/or  bomb  weapon 
deliveries.  In  gunnery,  it  is  the  harmonization  point;  i.e.,  where  the  gun 
will  shoot  at  harmonization  range,  usually  2250  feet  in  front  of  the  parked/ 
jacked-up  aircraft.  The  gunnery  ADL  is  usually  identifiable  by  either  an 
etched  gun-cross  reticle  or  the  caging  (Zero  Sight  Line)  of  the  movable- 
reticle.  Rarely  is  the  zero-sight  and/or  ADL  displaced  from  the  body  sym- 
metrical plane.  Thus,  the  FSRF/ADL  RFs  are  generally  small  or  no  y-axis 
rotations  of  the  BRF. 

The  lead  computing  gyro  case/caged  RF  (LCGCRF)  defines  zero  gunsight 
deflection.  While  some  older  sights  rotated  the  case  about  y for  rolling 
sensitivity,  the  x-axis  of  some  of  the  newer  LCGCs  is  aligned  parallel  to 
the  gunnery  ADL. 

The  gun  RF  is  similarly  defined  except  it  is  rotated  about  both  y and 
z axes  of  the  BRF  in  harmonization  to  the  ADL  and  also  to  correct  for  both 
gravity  drop  and  parallax.  Note,  no  roll  is  defined,  irrespective  of  gun 
receiver  block  orientation.  One  exception  to  parallax-gravity  correction 
is  noted  where  the  gun  fires  parallel  to  the  BRF  x-axis  and  all  corrections 
are  made  within  the  gunsight  mechanization.  The  ADL  is  parallel  to  the  FRL 
in  this  instance. 

It  is  recommended  that  all  of  the  above  definitions  be  in  the  dash  V) 
Technical  Orders  for  USAF  aircraft  of  interest. 


SECTION  V 


AIRCRAFT  VARIABLE  REFERENCE  FRAMES 

The  principal-axes  RF  (PRF)  defines  a set  of  axes  for  which  guidance 
control  inputs  most  nearly  cause  couplets,  pure  rotations  of  the  airframe. 
This  simplifies  flight  equations  by  uncoupling  the  axes;  i.e.,  rotation 
about  one  axis  does  not  affect  rotations  in  the  other  two.  In  aerodynamic 
terminology,  the  cross-product  inertias  (CPIs)  are  zeros.  It  is  in  this  Rl 
that  the  primary  aerodynamic  coefficients/stability-derivatives  are  defined. 
01  course,  the  PRF  varies  with  flight  condition  and  configuration.  Most 
modeling  ignores  this  RF , because  flight  is  also  defined  for  the  stability 
axes  Rl  (SRI);  i.e.,  the  coefficients  being  derived  for  the  SRF. 

The  SRF  is  rotated  from  the  body  by  the  wind  angle  components,  alpha 
(AOA)  and  beta.  The  SRF  is  identical  to  the  wind-axes  RF  (WRF)  for  steady- 
state  flight  conditions  and  is  the  basis  of  development  of  the  smal 1 -pertur- 
bation theory-derived  disturbed-axis  RF  (DRF) . The  DRF  is  identifiable  as 
the  time  varying  WRF/SRF  output  of  the  model;  transformations  being  made 
with  respect  to  the  instantaneous  values  of  the  WRF/SRF.  It  is  noted  that 
there  are  beta  (negative  psi)  and  alpha  (theta)  rotations  between  the  BRF 
and  WRF/SRF  but  no  roll.  Note,  CPIs  arc  non-zeros,  but  most  models  ignore 
them.  Of  course,  good  flight/stability  control  augmentation  systems  are 
designed  to  compensate  for  CPIs.  Accordingly,  omissions  are  tolerable  within 
1 imitations . 

Similar  to  the  FSRF,  the  dynamic  sight  pipper  and  target  line  of  sight 
(LOS)  RF  (PLRF  and  TLRF)  originate  on  the  sight  combining  glass  and  like 
the  WRF/SRF  have  no  roll.  They  are  expressly  defined  by  the  sight  elevation 
(pitch)  and  azimuth  (psi)  deflections  of  their  respective  LOS  with  respect 
to  the  zero  sight  line.  Pipper  and  target -centered  impact  plane  definitions 
stem  from  these  RFs. 


SECTION  VI 

EXTERIOR  BALLISTICS  REFERENCE  FRAMES 

lhe  Applied  Research  Laboratories  (ARL)  of  the  University  of  Texas 
authored  References  3 and  4,  using  unique  indexing  of  the  RI-'s  which  requires 
mixed  sensings.  Section  II  indicated  that  this  causes  unnecessarily  compli- 
cated bookkeeping.  Also,  the  same  RF  is  denoted  differently  in  the  two 
references  and  in  various  sections  of  Reference  3.  The  most  distinctive 
notations  were  selected  for  presentation  herein. 

Figure  I shows  the  differences  in  indexing  of  the  basic  and  ARL  RFs. 

The  ARL  inertial  RF,  denoted  XI,  ETA,  ZliTA  ( / , n,  V)  , is  rolled  270  degrees 
with  respect  to  the  basic  RF , apparently  done  for  the  sake  of  calling  up 
positively,  initialized  with  ETA  directed  opposite  to  the  gravity  vector 
and  XI  directed  along  the  bullet  launch  heading.  Heading/elevation  about 
the  ARL  Y/Z  axes  is  oppositely/directly  sensed  to  that  of  the  basic  RF.  In 
the  figure,  N/E/FWD/RT/DWN  have  been  used  to  abbreviate  North/East/forward/ 
right /down. 

It  is  evident  that  the  ARL  author  is  consistent  with  indexing  the  RFs. 
The  reader  is  cautioned  that  in  following  discussions,  the  equivalent  ARL 
RF  has  the  270-degree  roll  difference  in  indexing  from  the  basic  which  will 
not  be  explicity  noted  for  each  RF. 

The  trajectory  initialization  RF  (TIRF),  denoted  1',  J',  k',  by  ARL  is 
defined  as  the  translation  (PPL1 , projectile  launch  position  vector  in  the 
FIRF)  of  the  gun  muzzle  at  fire  and  launch  heading  (ip- sub  L)  and  elevation 
(9-sub  L)  rotations  to  initial  bullet  direction  in  the  FIRF.  The  TIRF 
establishes  the  basis  for  which  a selection  from  the  various  ballistics 
models  can  be  used  to  model  the  trajectory. 

It  is  noted  that  TIRF  X-FWU  is  identically  the  P axis  of  the  Modified 
Point  Mass  (MPM)  P-Q  RF  which  is,  of  course,  the  RF  for  which  ARL  develops 
both  the  iterative  MPM  and  Closed-form  Siacci  Solutions  (CSS);  the  MPM 
being  the  apparent  USAF  ATAGAM  standard  although  ARL  suggests  that  the  CSS 
is  the  ideal.  Q is  GVD  from  the  tip  of  the  P vector  for  clarity  in  vis- 
ualizing bullet  location,  but  the  vector  analyst  can  define  Q from  the 
TIRF  origin  with  no  loss  of  generality.  While  the  MPM/Siacci  RF  is  single/ 
vertical  plane  only,  swerve  jump  components  can  be  normal  to  this  plane. 

This  is  sensed  for  the  TIRF  y-right  definition  for  consistency. 

The  following  definitions  for  b degrees  of  freedom  (6P0F)  models  are 
included  for  the  modeler's  general  understanding  and  appreciation  of  the 
approx i~>t ions  involved,  although  it  will  be  made  apparent  that  they  are 
not  used  directly  in  aerial  gunnery  modeling.  It  is  noted  that  the  ARL 
author  refers  to  "approximate  equations  of  motion"  for  6D0F  models  which 
is  true,  but  without  explaining  that  they  are  the  most  accurate. 
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Flying  out  the  trajectory,  the  6D0F  model  requires  a trajectory  flight 
RF  (TFRF) , denoted  1,  J,  K by  ARL,  analogous  to  the  wind/stability  (axes) 

RF  (WRF/SRF)  of  Section  V.  It  is  indexed  with  respect  to  the  TIRF  and  up- 
dates/det incs  relative  bullet  eg  translations  and  rotations. 

The  bullet  body  (BB)  RF,  defined  A,  B,  3 by  ARL,  is  essentially  identi- 
cal to  the  TFRF  excepting  BB  roll  rotation  after  the  yaw  on  the  bullet  damps 
out.  Thus,  the  effect  of  yaw  is  the  essential  differentiation  of  the  6D0F 
from  simplified  models.  Note,  the  BBRF  Y/ARL-3  and  Z/ARL-B  rotate  at  the 
BB  spin  rate  about  X/ARL-A.  The  BBRF  is  initialized  to  the  TFRF/TIRF  by  the 
rotations,  0-sub  Y and  ip-sub  Y which  are  defined  in  the  next  section.  A 
yawed  BBRF  X/ARL-A  precesses/nutates  about  the  TFRF  X/ARL-1.  To  account  for 
its  orientation,  ARL  defines  a 1,  2,  3 RF,  a bullet  nutation-accounting  RF 
(BNRF) , where  1 is  velocity/flight  vector  directed;  i.e.,  1 is  the  same  in 
the  1,  2,  3 and  1,  J,  K RFs.  Then  nutation  of  2 from  J (or  3 from  K)  defines 
the  precession  angle,  ARL  <p,  while  the  magnitude  of  the  yaw/sigma  angle, 

ARL  6,  is  defined  from  1 to  A in  the  1,  2 plane. 

It  is  noted  that  the  MPM/CSS  requires  only  the  definition  of  the  TIRF 
in  the  ATAGAM  FIRF.  The  6D0F  model  also  requires  the  TFRF,  BBRF,  and  BNRF. 
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SECTION  VII 
TRANSFORMATIONS 


Sections  4 and  5 of  Reference  2 present  a clear,  concise  treatment  of 
the  standard  transformations.  However,  methods  of  uncoupling  the  various 
parts  may  not  be  readily  evident  to  the  unfamiliar.  For  example,  the  I*  of 
subsection  5.1  could  be  related  to  projectile/bullet  position,  the  0 to  the 
body  eg,  and  the  01  to  the  origin  of  the  F1RF.  While  mathematically  correct, 
the  development  would  readily  become  mind  boggling  in  definition  of  the  X, 

Y,  and  2 derivatives  of  P with  respect  to  0. 


in  practice,  each  aircraft  and  bullet  eg  is  initiated  and  updated  in 
the  F1RF,  because  the  equations  of  motion  are  both  independent  and  defin- 
able in  this  RF.  For  each  axis,  both  translational  and  rotational  motions 
are  updated  by  the  Taylor  Series  equation: 


f(t)  = f(o)  + f • t + f * 1 2/2 ! + f • t 3 / 3 ! + HOTs 


(I) 


where  f = x/y/z  translation  and  <J)/0/i|)  rotation 

t = time  from  initiation,  zero  time 
HOTs  = higher  order  terms  (usually  neglected). 

As  previously  noted,  the  bullet  equations  of  motion  (excluding  jump)  are 
usually  defined  in  a single  plane  and  then  transformed  to  the  FIRF. 

It  is  noted  that  the  time  increment  in  iterative  models  is  kept  small 
so  that  the  HOTs  have  negligible  effect.  Many  ignore  the  third-derivative/ 
jerk  term.  Accuracy  is  affected  by  both  iteration  interval  and  truncations 
of  the  series.  The  resultant  angles  are  Euler  attitudes  of  the  body  RF 
with  respect  to  the  MIRF  indices. 

In  the  real  world,  the  Euler  attitudes  are  sensed  by  pickoffs  on  the 
gyro  gimbals.  If  inertial,  the  platform  contains  rate  gyros  and  (trans- 
lational) accelerometers.  However,  the  rate  gyro  and  accelerometer  outputs 
are  used  internally.  That  is,  the  normal  outputs  of  the  platform  are 
attitudes  and  translations.  Thus,  models  that  use  real  world  inputs  do  not 
require  accountings  of  time  and  derivatives  to  determine  Euler  attitudes. 
However,  instrumentation  is  most  often  found  lacking  and  occasions  arise 
where  angular  derivatives  in  the  body  RF  must  be  obtained.  Of  course, 
differentiation  is  required  with  the  well  known  noise  problems.  Most 
models  define  body  rates  and  Euler  attitudes  must  be  obtained  from  them. 
Thus,  the  required  rotational  transformations  (RTs)  can  be  expressed  as: 
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where  S,  l',  T prefixes  are  shortening  abbreviat ions  for  sine,  cosine,  and 
tangent  functions.  Due  to  sequential  rotations  in  their  definition  (Kef 
ercnce  2),  these  RTs  are  not  orthogonal.  Properties  of  orthogonality  are 
defined  later  in  this  section. 


Before  general  use  of  large  digital  processing  machines,  limited/analog 
processors  required  linearizing  p/q/r  to  phi-dot/theta-dot/psi-dot  for 
specified  initial  conditions;  i.e.,  near  zero  psi,  theta,  and  phi  where 
sine/cosine  functions  are  approximately  0/1.  The  restrictions  on  lineariza- 
tions are  obvious. 


Because  the  Hu  1 or- sensing  gyro/platform  is  affixed  to  the  body,  occasions 
arise  quite  frequently  where  understanding  that  filler  and  body  rates  arc  not 
identical  is  difficult  for  some.  Note  that  the  M1RF  axes  remain  inert  ially 
fixed  and  sense  any  of  the  arbitrarily  oriented  body  motions  in  their  re- 
spective axes. 


Having  established  the  relationships  between  filler  and  body  axes  rates 
and  the  definition  of  the  filler  attitudes,  the  standard  RTs  between  RTs  can 
be  expressed  as: 
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whore  St,C  prefixes  are  shortening  abbreviations  for  sine  and  cosine  functions, 
RX/RY/RZ  is  the  RT  about  the  X/Y/Z  axis,  A/B/C  is  the  associated  Rll  rotation 
angle  from  the  reference  index.  It  is  noted  that  the  development  allows  in- 
dependent determination  of  components  in  each  rotation  plane  (i.e.,  the  XY 
plane  for  Z-axis  rotations,  etc.)  separately,  while  the  results  for  all  of 
the  ordered  rotations  are  given  by  the  ordered  product  of  their  RTs. 


The  majority  of  the  objects  of  interest  requiring  transformation  are 
associated  with  the  attacker  airframe  whose  states  must  be  updated  to  define 
flight,  while  the  only  object  of  interest  in  the  WRF  is  the  velocity  vector. 
Thus,  while  object-to-body-to-wind-to-earth  and  reverse  transformations  can 
be  defined  correctly,  they  contain  needless  complications.  Relating  objects 
to  the  body  and  using  body-to-earth  and  reverse  transformations  minimizes 
the  number  of  transformations  required.  Accordingly,  the  standard  trans- 
formations are  presented  as  follows. 

i.ARTII- TO- BODY  AND  REVERSE 

PQB  = [RX(4>)  • RY  (0)  • RZ(i(0]  • PQE  = R BE  • PQE  (7) 

where  PQB/PQE  is  the  body/earth  column  matrix  of  X,  Y,  and  Z position 
coordinates/components  of  some  point  of  interest  (Q)  such  as  the  target, 
projectile,  gun  muzzle,  sight,  etc.  The  psi,  theta,  phi  ordering  in  matrix 
algebra  is  evidenced  in  the  pre-multiplications. 

P/V/A  first-letter  notation  identifies  position/velocity/ acceleration 
(.vectors).  The  intermediate  letter(s)  identify  the  point/object  of  interest, 
while  the  last  letter  identifies  the  RE.  Coordinates/components  are  iden- 
tifiable by  the  X,  Y,  or  Z ending  added  to  the  vector  notations.  Elements 
ot  the  RTs  are  conventionally  identified  by  first  letter  C,  second  letter 
symbol  for  object,  third  letter  notation  for  output  RE  axis  component,  and 
fourth  letter  for  input  RE  axis  contributing  component. 

Using  the  orthogonality  principle  of  RTs, 


RBA 1 = UBA  1 = RAB 


(8) 


where  \ is  one  RE,  B is  another,  and  AB  implies  transforming  from  B to  A 
coord i nates , 


PQE  RBI  1 • PQB  - R!  B • PQB 


(9) 
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The  AB  notation  used  above  is  totally  reversed  or  used  inconsistently  by 
some  modelers.  Consistency  with  that  notation  has  been  used  in  the  following 
discussions.  It  is  also  noted  that 

HI  (II)  * = Rl(-D)  where  D = associated  angle  for  any/1  axis  (10) 


and 

RBL1  = RZ(-ip)  • RY  (-0)  • RZ(-4>)  (11) 

i.e.,  rotating  backwards  through  the  angular  transformations  is  the  same 
as  the  transpose. 

Standard  notation,  proffered  by  and  received  from  AFAL/RWT-2  in  coor- 
dination, is  repeated  verbatim  as  Appendix  A of  this  report.  Note  reverse 
AB  implication.  Coordinated  agreement/refinement  and  a listing  of  common 
variables  is  considered  to  be  a good  basis  for  a companion  report. 

WIND -TO- BODY  AND  REVERSE 

Flight  definitions  define  the  wind  vector  and  angular  components  alpha 
and  beta  from  the  wind  to  either  the  body  or  other  FBRF.  Usually,  the  angles 
between  the  WCRF/other  FBRF  frame  and  the  body  are  available  so  rather  than 
using  multiple  transformations,  the  angles  to  the  body  can  be  obtained  by  a 
simple,  directed-sum.  Thus,  the  basic  wind/body  transformations  arc: 

PQB  = [RY (a)  • RZ(-S) ] • PQW  = RBW  • PQW  (12) 

and 

PQW  = RBW1  • PQB  = RWB  • PQB  (13) 

where  PQW  is  defined  similarly  to  PQB/PQE/  for  X,  Y,  2 components  in  the  wind 
axes.  Note  that  accounting  for  left-hand  sensed  beta  is  accomplished  by 
inputting  negative  beta  to  the  RHRT.  Note  also  the  absence/identity  of  the 
roll  (x-axis)  RT. 

EXTERIOR  BALLISTICS 

There  arc  several  factors  that  are  required  to  initialize  the  fired 
bullet;  i.e.,  the  launch  velocity,  its  position  and  orientation,  and  the 
wind  angle  on  the  projectile  if  windage  jump  is  considered. 

The  essential  launch  velocity  vector  (VI,)  transforms  are  wi nd-to-body 
aircraft  velocity  vector  (VA) , addition  of  muzzle  velocity  vector  (VM)  and 
muzzle  tangential  velocity  vector  (VMT)  in  the  same  frame,  and  rotating 


their  sum  into  the  FIRF.  VA  in  the  BRF  can  be  expressed  as: 


VAB  = RBW  • VAW  (14) 

where  VAB/VAW  is  a column  vector  of  the  X,  Y,  and  Z velocities  in  the  BRF/WRF. 
The  components  in  the  WRF  are  VA,  0 and  0;  so  that  RBW’VAW  can  be  reduced 
to  the  first  column  times  VA  (scalar)  if  desired  for  computational  efficiency. 

The  muzzle  velocity  vector  in  the  BRF  can  be  expressed  as: 

VMB  - IRY(FC)  • RZ (AG) ] 1 • VMG  = RBG  • VMG  (15) 

where  l.G/AG  is  the  theta/psi  rotations  of  the  gun  from  the  BRF.  Again,  VMG 
components  in  the  gun  RF  are  VM,  0,  0;  so  RBG*VMG  can  be  reduced  if  desired. 

The  gun  muzzle  and  sight/LCG  are  usually  located  at  significant  distances 
from  the  aircraft  eg.  Aerodynamic  flexures  are  generally  unknown  quantities 
so  that  these  elements  are  considered  fixed  in  the  airframe.  Since  a point 
on  a rigid  arm  rotates  with  the  center/cg,  the  sight/LCG  offsets  do  not 
affect  their  functionings.  However,  the  tangential  velocity  of  the  muzzle 
adds  vectorially  to  the  launch  velocity.  The  tangent  velocity  vector  in  the 
BRF  can  be  expressed  as: 

VMTB  = w@R  = 2 • PMB  (16) 

where  PMB  is  the  muzzle  position  vector  in  the  BRF,  and 
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The  tilde  operator  allows  cross-product  operations  to  be  defined  by  matrix 
multiplication.  Thus,  the  components  are: 


VMTBX  = qz 

- ry 

(17) 

VMTBY  = rx 

- P- 

(18) 

VMTBZ  = py 

- qx 

(19) 

To  ballpark  the  effect  of  tangential  velocities  on  highly  maneuverable  air- 
craft, the  upper  limits  of  p/q/r  are  about  4/0.55/0.1  rad/scc  while  typical 
muzzle  displacements  are  10-20/0- 10/0-5  feet.  Accordingly,  comparisons  of 
tangential  velocities  of  a few  ft/sec  to  average  bullet  velocities  of  a few 
thousands  ft/sec  indicate  an  effect  of  about  a thousandth  of  a rad/sec  change 
in  direct ion/t ime-of-fl ight  (TOF) . Because  of  their  small  effect,  some 
models  ignore  tangential  velocities. 
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The  total  launch  velocity  in  the  BRF  and  F1RF  can  be  expressed  as: 


VLB  = VAB  + VMB  + VMTB 


VLI  = REB  • VLB  (21) 

It  is  noted  that  the  vectorial  summing  in  these  equations  inherently  accounts 
for  what  bal listicians  have  called  aerodynamic  jump;  i.e.,  the  vector  sum  of 
aircraft  and  muzzle  velocity  vectors. 

The  projectile  position  at  launch  in  the  FIRF  can  be  expressed  as: 

PPLI  = PA I + REB  • PMB  (22) 

where  PAI  is  the  inertial  firer-cg  location  vector,  defined  by  Equation  (1). 

The  wind/yaw  angle  on  the  projectile  is  simply  the  directional  difference 
of  the  launch  velocity  vector  from  the  muzzle  velocity  vector.  The  muzzle 
velocity  vector  in  the  FIRF  is  simply: 

VMI  = REB  • VMB  (23) 

Thus,  the  launch  and  muzzle  velocity  vector  directions  can  be  expressed  as: 

ifi-sub  L = arctan  (VLIY/VLIX),  (24) 

iji-sub  M = arctan  (VMIY/VM1X) , (25) 

6-sub  L = arcsin  (-VLIZ/abs(VLI) ) , (26) 


0-sub  M = arcsin  (-VMIZ/abs (VMI) ) , (27) 

where  abs  (A)  = magnitude  of  A,  i.e.,  the  root-sum-square  of  the  components. 
Accordingly,  the  bullet  yaw  components  can  be  expressed  as : 

4)-sub  Y = (ip-sub  M)  - (tp-sub  L)  (28) 


6-sub  Y = (0-sub  M)  - (6-sub  L) 

The  spatial  bullet  yaw  (E)  can  be  defined  by  either: 

abs (VLI) abs (VMI)  cos  (E)  = VLI Q VMI 


m 


or 

E = arc  cos  (cos(ip-sub  Y)  cos  (0-sub  Y))  (31) 

The  high  dynamics  of  precessional/nutational/corkscrewing  motion  of  the 
yawed  projectile  requires  an  extremely  small  iteration  interval,  even  though 
the  ballistics  flight  equations  are  very  simple  in  comparison  to  aircraft. 
Thus,  even  CPC  (>(>00  processor  time  is  about  100  times  real-time.  Because 
the  corkscrew  essentially  damps  out  in  the  first  500  feet  of  the  aerial 
trajectories  of  interest  and  its  disturbance  effects  arc  predictable  from 
analysis  of  <>IK)I  outputs,  use  of  such  modeling  in  general  analyses  would  be 
i ne  ft"  i c i ent . 

The  effects  of  yaw  are  obviously  direct  lift,  niagnus  lift,  and  increased 
drag.  ARL  has  shown  in  References  3 and  4 that  increased  drag  can  be  mech- 
anized by  the  addition  of  an  exponential  term  in  the  velocity  equations. 
Conversion  to  basic  RFs  requires  consideration  of  different  indexing  in- 
cluding the  delta  altitude  term.  Thus,  the  remainder  of  yaw  effects  defi- 
nition involves  windage/swerve  jump. 

Figure  5 shows  the  geometry  of  windage/swerve  jump  with  the  TFRF  X/ 
velocity-vector  into  the  page.  The  yaw/sigma  (E)  angle  from  the  flight 
vector  to  the  projectile  nose  direction  is  represented  by  the  radius  of  the 
circle.  The  resultant  effects  over  the  corkscrew  path  of  the  projectile 
can  be  expressed  as  a radially  outward  effect  (vector  a)  and  a tangential 
effect  (vector  b)  with  regard  to  the  projectile's  initial  orientation.  If 
the  radial  and  tangential  effects  were  directly  correlatable  with  the  initial 
orientation,  then,  the  rad i a lly-out ward/ tangent ia 1 effect  could  be  attrib- 
uted to  direct/magnus  lift.  The  drawing  of  the  vectors  on  top  of  the  circle 
arc  included  for  ease  in  visualizing  the  cause  and  effect,  while  their 
orientation  on  the  lower  right  of  the  circle  are  included  to  illustrate 
arbitrary  orientation  of  initial  conditions,  defined  by  ARL  <J>.  The  mag- 
nitudes of  a and  b are  derived  as  a function  of  P (which  is  closely 
correlatable  to  projectile  range)  for  specified  initial  orientations.  Thus, 
even  if  there  is  a phasing  between  initial  orientation  and  effect,  the  values 
of  a and  b will  reflect  the  effect.  Thus,  it  follows  directly  that: 

delta-y  = P • E (a  sin  4>  + b cos  <}>)  (32) 

and 

delta-z  = -P  • E (a  cos  <J>  - b sin  i|>)  (33) 

where  P = Siacci  space  function. 


It  is  noted  that: 


' 


and 


so  that 


and 


E sin  <j>  = i|/-sub  Y 


delta-y  = P (a  • tp-sub  Y + b • 0-sub  Y) 


(35) 


(36) 


dclta-z  = -P  (a  • 0-sub  Y - b • ip-sub  Y)  (37) 

The  modified  Alii,  estimates  of  a and  b are  0.0  and  0.037  rad/ rad  for  M-50 
series  20  nun  projectiles.  Thus,  the  terms  involving  a in  Equations  (36) 
and  (37)  could  be  dropped  for  this  special  case.  The  general  forms  of 
these  equations  are  desirable  for  compatibility  with  other  projectiles. 


Referring  to  Figure  6,  it  can  be  seen  that  an  ambiguity  arises  in  trans- 
forming delta-Z  into  the  P and  Q axes.  Obviously,  delta-Z  times  the  cosine 
of  0-sub  L is  the  contribution  to  Q which  is  GVD.  But,  by  definition 
delta-Z  is  normal  to  P.  Specifically,  at  9-sub  L of  90  degrees,  no  component 
of  delta-Z  can  be  attributed  to  either  Q or  P.  Thus,  the  first  step  is  to 
transform  P-Q  coordinates  to  TIRF  coordinates.  The  T1RF  coordinates  can  be 
expressed  as: 

X = P - Q sin  (0-sub  L)  (38) 


Y = del ta-Y  (39) 

Z = Q cos  (0-sub  L)  + delta-Z  (40) 

The  transformation  from  the  TIRF  position  vector  (PPT)  into  the  FIRF 
vector  (PPI)  can  then  be  expressed  simply  as: 


PPI  = PPLI  + REB  • PPT 


(41) 


where 


REB  = (RY  (0-sub  L)  • RZ  (i|>-sub  L))^. 

It  is  noted  that  Equation  (41)  establishes  the  position  but  not  the 
direction  of  the  bullet  in  the  FIRF.  Precise  bullet  orientation  is  impor- 
tant only  to  precise  impact  definition  which  is  generally  not  defined;  i.e., 
target  surfaces  are  not  specifically  defined.  Accordingly,  most  models 
ignore  the  gravity  rotation  which  is  generally  less  than  a degree  and 
windage  jump  components  which  are  generally  about  a few  mi  1 1 iradians . Deri- 
vation of  bullet  direction  is  relatively  simple.  The  TFRF/TIRF  Y is  always 
in  a level  FIRF  plane  so  that 
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Horizontal 


u)-sub  Y = -V-sub  B/(g  • cos  (0-sub  B) 


(42) 


where  ca-sub  Y is  the  gravitational  rotation  rate 


V-sub  B = bullet  velocity  in  ft/sec 


2 

g = gravitational  acceleration  ft/sec 

Starting  at  launch  where  0-sub  B is  equal  to  0-sub  L,  0-sub  B can  be  updated 
by  the  Taylor  Series  equation.  Finally,  the  jump  angle  components, 

0-sub  .J  = -delta-Z/P  and  i|i-sub  J = delta-Y/P,  are  directly  additive  to 
U-sub  B and  ijj-sub  B ( = i|>-sub  L)  , respectively,  as  a first  order  approx- 
imation, i.e.,  the  change  in  gravitational  rotation  due  to  0-sub  J is  a 
second-order  effect. 

It  has  been  made  apparent  that  6D0F  models  are  inefficient  in  general 
ATAGAM,  due  to  excessive  running  time.  In  iterative  solutions,  accuracy  is 
a function  of  iteration  interval  which  should  be  apparent  from  basic  calculus. 
For  general  aerodynamic  flight,  a 0.1-second  interval  is  considered  to  be 
the  best  compromise  between  accuracy  and  running  time.  The  approximation 
is  evident. 

BULLET-TO-SIGHT 

Quite  a few  models  compare  bullet,  target,  and  pipper  coordinates  in 
the  FSRF.  Since  the  target-to-sight  transforms  are  of  identical  form,  only 
the  bul let-to-sight  transforms  are  presented. 

It  is  noted  that  the  aircraft  eg  is  in  motion  with  respect  to  the  bullet 
launch  position  and  coordinates  in  the  FSRF  are  relative  ones.  Thus,  the 
first  transformation  is  to  derive  relative  coordinates  of  projectile  position. 
The  bullet-position-sight-relative  vector  (PPSI)  of  inertial  components  with 
respect  to  the  FSRF  can  be  expressed  as: 

PPSI  = PPI  - PA I + REB  • PSB  (43) 

where  PSB  is  the  FSRF-origin  location  vector  in  the  BRF. 

The  transformation  to  the  BRF  can  be  expressed  as: 

PPSB  = RBE  • PPSI  ; (44) 

and,  finally,  the  transformation  to  the  FSRF  can  be  expressed  as: 

PPS  = RSB  • PPSB  (45) 
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where  RSB  is  RY  (0-sub  IS);  0-sub  FS  being  the  rotation  of  the  FSRF  from 
the  BRF.  Of  course,  abs  (PPS)  is  the  relative  projectile  range,  and  the 
angular  positions  with  respect  to  the  FSRF  origin  in  the  FSRF  can  be 
expressed  as: 


i^-sub  P = arcsin 

(PPSY/abs  (PPS) ) 

(46) 

0-sub  P = arcsin 

(-PPSZ/abs  (PPS) ) 

(47) 

Presentation  of  all  projectile  coordinates  at  a specific  instant/iteration 
is  analogous  to  the  visibly-sensible  tracer  line.  Sequential  presentation 
of  their  intersection  (requiring  interpolation  of  fly-by)  in  the  impact 
plane  is  analogous  to  what  would  be  seen  if  the  projectiles  fuzed,  regardless 
of  hit,  at  fly-by.  The  latter  is  more  important  in  modeling  hit  probabilities. 

The  target,  projectile(s)  (at  impact-plane  intersection),  and  the 
pipper  positions  in  the  sighting  RF  allow  calculations  of  accuracies  and 
statistics. 

ANGLE -OFF -TRANSFORMAT IONS 

In  modeling,  angle-off  (AO)  is  defined  to  be  the  angle  between  the  LOS 
range  vector  (LOSR)  and  the  target  velocity  vector.  In  the  past,  AO  has 
been  defined  as  an  acute  angle,  due  to  geometric  limitations  on  accurate 
gunsight  operations.  The  following  derivations  have  been  generalized. 

The  position  of  the  attacker's  sight  (PSI)  in  the  F1RF  is  simply  the 
last  two  terms  of  Equation  (43) ; 

PSI  = PA I + REB  • PSB  (48) 

Then,  the  FIRF  LOSR  (PLI)  can  be  expressed  as 

PLI  = PT1  - PSI  (49) 

where  PTI  is  the  target-cg  location  in  the  FIRF.  It  is  noted  that  angle-off 
is  defined  for  the  conical  angle  of  PLI  about  the  target  velocity  vector; 
i.e.,  any  AO  that  is  defined  by  a point  on  the  cone  is  not  distinguished 
from  another.  Thus,  AO  can  be  defined  by  dot  producting  PLI  with  the  tar- 
get velocity  vector  (VTI)  in  the  FIRF;  i.e., 

cos  (AO)  = PLlOvTI/(ABS(PLI)  • abs  (VTI))  (50) 

It  is  also  noted  that  angle-off  is  usually  considered  a cursory  measure  and 
most  models  neglect  the  PSB  term. 
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Another  term,  angle-off-the-tai  1 (AOT)  is  used  as  a gross  measure  of 
target  aspect.  AOT  is  sometimes  defined  as  the  negative  of  PLI  dot-producted 
with  a unit  vector  (UV)  along  the  target  negative  X axis  which,  of  course, 
is  identical  to  the  dot  product  of  the  position  vectors;  i.e., 

UVTXI  = RFB  • (1,0, 0)F  (51) 

where  UVTXI  is  the  FIRF  transform  of  the  TBRF  (target  BRF)  X UV. 

1 is  the  X component  on  the  UV  along  the  TBRF  X 


and, 


cos  (AOT)  = (PLI O UVTXI) /abs (PLI)  (52) 

Note  that  RLB  is  defined  for  the  TBRF  in  this  equation. 

In  some  models,  target  aspects/views  are  defined,  using  the  basic  psi, 
theta,  phi  conventions.  This  requires  definition  of  the  negative  (PNLI)  of 
PLI  in  the  TBRF  and  then  defining  the  directions  of  the  PNLI  in  the  TBRF. 
This  can  be  expressed  as: 


PNLB  = RBE  • PNLI  (53) 

where  PNLB  is  the  PNLI  vector  in  the  TBRF.  Again,  RBF.  is  defined  for  the 
TBRF. 

Then,  the  target -aspect  heading  (HASP),  pitch  (PASP) , and  roll  (RASP) 
angles  can  be  expressed  as: 


HASP  = arctan  (PNLBY/PNLBX)  (54) 
PASP  = FCN1  (HASP,  PARG)  (55) 
PARC  = arcsin  (PNLBZ/abs (PNLB) ) (56) 
RASP  = arctan  (PNLBZ/PNLBY)  (57) 


where 


PARG  is  the  pitch  aspect  angle  argument 

FCN1  = PARG  for  HASP  in  forward  hemisphere 
180-PARG  for  HASP  in  rear  hemisphere. 

The  ambiguities  for  HASP/RASP  at  ± tt/2  can  be  handled  in  any  consistent 
manner,  being  usually  assigned  to  the  positive  semi-circle.  It  is  noted 
that  aspects  in  the  vertical  plane  of  the  target  define  positive  90-degree 
roll  aspects  which  should  not  be  confused  with  target  attitude  in  the  FIRF. 


Also,  note  that  rotation  of  pilot/firer  up  about  the  LQSR  does  not  change 
these  target  aspects.  That  is,  whether  the  firer's  head  is  upright,  tilted, 
or  inverted  in  the  FIRF  does  not  change  the  fact  that  he  views  the  target's 
top/bottom,  left/right  side,  and  front/rear  views/aspects  from  the  specified 
aspect  angles. 

Finally,  crossing  angle  is  a term  that  is  used  to  define  the  bullet 
aspect  with  respect  to  the  target  fly-by.  It  is  identically  derived  to 
AOT  and/or  aspects  with  the  exception  of  using  the  negative  of  the  bullet 
velocity  vector  instead  of  PNLI.  Of  course,  striking  angles  against  and 
unit  projected  areas  of  front,  right  side,  and  bottom  views  are  defined  by 
dot-producting  the  negative  of  the  bullet  velocity  vector  (VNB)  with  UVs 
along  the  TBRF  X,  Y,  and  Z axes,  respectively.  Note  that  negative  products 
of  the  UV  and  the  VNB  component  in  its  direction  indicate  rear,  left-side, 
and  top  views.  The  total  projected  area  of  a view  is  simply  the  normal  area 
of  the  view  times  its  unit  projected  area.  Where  vulnerabilities  of  the 
target  are  given  by  azimuth  and  elevation  aspects  with  elevation  as  the 
combination  of  roll  and  pitch  aspects,  the  complement  of  the  striking  angle 
against  the  top/bottom  view  can  be  used  to  define  this  elevation. 

MAJOR  CONFUSION  FACTORS  IN  TRANSFORMATION 

Confusion  rarely  occurs  when  it  is  made  clear  that  coordinates  of  a 
specific  object  are  being  transformed  from  one  coordinate  system  to  another. 
In  many  cases,  the  model  develops  multiple  transformations  without  reference 
to  what  is  to  be  transformed  and  to  which  RFs.  Also,  some  models  will  refer 
to  Euler  attitudes  as  Earth  angles  rather  than  body  angles  to  Earth  indices/ 
references.  Specific  comments  of  what  is  being  transformed  nd  the  RTs 
involved  will  avoid  such  confusion. 


RirOMMIiNDATIONS  AND  CONCLUSIONS 


1'he  preceding  discussions  are  directed  to  provide  a common  basis  for 
development  of  ATAUAM.  Compliance  with  these  preferred  standards  should 
simplify  their  structure  and  general  comprehension. 

Ihe  discussions  understandably  fall  short  of  total  comprehension.  A 
first  addendum  may  be  envisioned  to  expound  on  the  various  models  used  to 
define  flight.  This  is  the  major  difference  between  most  tactical  models. 

The  simpler  models  consider  flight  control  only  to  rotate  the  flight  vector, 
neglecting  direct  lift  effects  along  with  cross-coupling,  phugoid,  and  short 
period  dynamic  modes  of  the  airframe.  Another  could  address  the  paper  pilot, 
the  most  controversial  aspect  of  all-digital  codes,  besides  neuro-muscular 
lags  and  decision-delay  functions,  high  integral  gain  along  with  optimized 
proportional  and  derivative  gains  is  sometimes  included  in  the  definition 
of  a second-order  paper  pilot  to  mislead  the  unfamiliar  into  believing  that 
the  flight  model  encompasses  all  dynamics  of  t lie  airframe.  Others  could  be 
directed  to  explain  and  compare  various  interpolation,  integration,  differen- 
tiation, smoothing,  atmosphere,  matrix  algebra  routines,  etc.  that  make  up 
t he  standard  subroutines  of  the  model. 

It  is  intended  that  this  treatise  plus  other/addenda  will  provide  a 
basis  for  modular  construction  of  aerial  gunnery  models  whose  capabilities 
and  limitations  are  readily  identifiable. 


Rl  I HRl  NC'.l.S 
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APPENDIX  A 

RACEE  Simul at  ion -FORTRAN  Kinematic  Variables  Mnemonics 


The  FORTRAN  kinematic  variables  have  been  named  by  the  conventions 
below.  It  is  hoped  that  after  reading  this  appendix,  their  interpretation 
will  be  both  uncomplicated  and  precise.  The  intention  of  the  naming  con- 
vention is  to  distinguish  between  five  properties  of  these  variables.  These 
properties  are  as  follows: 

1.  Order  of  derivative  and  type  of  variable.  For  example,  position, 
velocity,  acceleration,  angular  velocity,  transformation  matrix,  etc. 

2.  Physical  object  or  mathematical  object  concerned.  For  example, 
attacker  aircraft,  target,  inertial  frame,  air  mass  frame,  line  of  sight 
frame,  etc. 

3.  Another  object  (physical  or  mathematical)  which  the  variable  of  the 
first  object  is  in  relation  to.  For  example,  velocity  (variable)  of  target 
(first  object)  with  respect  to  attacker  (second  object).  Generally,  the 
same  set  of  objects  could  be  used  here  as  in  Number  2. 

4.  The  coordinate  frame  into  which  components  of  the  vectors  are  re- 
solved. Noting  that  most  of  the  variables  in  Number  1 are  vectors,  and 
that  to  operate  numerically  on  a vector,  it  must  be  resolved  into  a coor- 
dinate frame,  these  frames  must  be  specified.  Hence,  for  a complete  de- 
scription of  the  variable  given  in  the  example  of  Number  3,  it  would  be 
necessary  to  specify  that  the  vector  was  resolved  into  attacker  body  axes, 
inertial  axes,  or  line  of  sight  axes,  etc. 

5.  It  is  recognized  that  there  are  other  attributes  that  might  be 
sytematized,  such  as  noise-corrupted  versus  perfect  variables.  This  could 
be  recognized  as  a fifth  property. 

The  mnemonics  are  determined  by  the  following  conventions: 


FORTRAN  name  = "uvwxyz" 


where  "uv"  = PS 
VL 
AC 
SF 
OM 


position 

velocity 

acceleration 

specific  force  = acceleration  - gravity 
omega  (angular  rate) 
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1 


01) 

omega  dot  (angular  acceleration) 

SI 

ip  (general  liuler  azimuth  angle  about  z-axis) 

Til 

0 (general  l:uler  elevation  angle  about  y-axis) 

I’ll 

4>  (general  liulcr  roll  angle  about  x-axis) 

T 

transformation  matrix. 

where 

w or  x or  y 

= A 

attacker  body 

T 

target  body 

I 

inertial  frame 

W 

wind  frame  of  attacker 

L 

line  of  sight  frame  to  target  (non-roll  stabilized) 

M 

missile  body 

B 

bu 1 1 et 

V 

bullet  velocity  frame  as  in  "MISS" 

S 

scoring  frame  as  in  "MISS", 

and  where 

z = 

(.blank)  uncorrupted  (exact)  program  variable 

M 

measured  varible,  e.g.,  true  variable  corrupted 
by  noise  for  input  to  airborne  algorithm 

Some  examples 

; of  this  convention  are: 

PSTAI 

target  distance  north,  east,  and  down  from  attacker 

VLAIA 

attacker  inertial  velocity  in  body  axes 

VbAWA 

attacker  air  mass  velocity 

0MA1A 

attacker  angular  rate,  i.e.,  p,  q,  r 

SITAL 

target  azimuth  angle  from  attacker 

T1A 

transformation  from  inertial  to  attacker  body  axes 

SFTIT 

target  specific  force  vector  (i.e.,  SFT1T(3) 

= -target  load  factor  • 32.174) 

S1A1A 

attacker  heading  angle  (of  body) 

One  of  the  benefits 
For  example,  in  any 

of  this  convention  is  rapid  checkout  of  vector  operations, 
addition  or  subtraction  of  vectors,  the  fifth  letter 
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must  be  the  same.  In  any  transformation  of  vectors,  the  fifth  letter  of  the 
transformed  variable  and  the  fifth  letter  of  the  original  vector  determine 
which  transformation  to  use.  An  example  of  this  would  be 

PSTAA  = (T1A)  (PSTA1 ) 

In  addition  or  subtraction,  the  third  and  fourth  letters  must  be  consistent 
with  the  operation.  For  example, 

VLTAI  = VLTI1  - VLAII 


or 


ACT 1 1 = ACTA I + ACAII 


IW II  , * 
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INITIAL  DISTRIBUTION 


Hq  USAF/SAMI  1 
Hq  USAFE/DOQ  1 
Hq  PACAF/DOOFQ  3 
Hq  TAC/DRA  1 
ASD/ENFEA  2 
ASD/ENESH  1 
AUL/LSE  71-249  1 
OO-ALC/MMWMP  2 
AFIS/INT  1 
DDC/DDA  2 
AFATL/DLODL  2 
AFATL/DL  1 
US  Army  TRADOC  Sys  Analysis  Act/ 

ATAA-SL  1 
ASD/XRP  1 
COMIPAC/ I -232  1 
Armament  Systems,  Inc  1 
AFAL/RWT  1 
ASD/XRT  1 
AMRL/MEB  1 
Hq  USAF/SAGF  1 
Hq  AFSC/SDZ  1 
NATC/SA53-A  1 
NSWC/DL  1 
NSWC/G-10  1 
USAMSAA/DRXSY -ADM  1 
USAMSAA/DRXSY - J 1 
6510  Test  Wg/TEEDM  3 
NWC/Code  318  1 
USAF  TFWC/TEM  1 
USAF  TFWC/OA  1 
Hq  SAC/DOOB  1 
Hq  SAC/NRS  1 
Center  of  War  Gaming  1 
Navy  Fighter  Weapon  School  1 
VX-4  1 
Fleet  Combat  Training  Ctr, 

Atlantic,  Code  221  1 
Surface  Warfare  Development  Gp/ 

Naval  Amphibious  Base  1 
Fleet  Combat  Direction  Service 

Training  Ctr,  Pacific  1 
AFATL/DLD  1 
AFATL/DLDG  1 
AFATL/DLY  1 
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AFATL/DLYV 

AFATL/DLYD 

ASD/SDX 

AFA/DFASC 

TAC/INAT 
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